Dynamic spin-glass behavior in a disorder-free, two-component model of quantum 

frustrated magnets 



(N 
O 

o 

(N 
S3 

a 
m 

(N 



CO 

c3 



S3 

o 

(N 
> 

O 
^1" 

Cs 
O 



c3 



S3 

o 
o 



X 
S3 



Frederic Mila (o) and David Dean^ 
Institut de Physique Theorique, Universite de Lausanne, CH-1015 Lausanne, Switzerland 
( b ' Laboratoire de Physique Quantique, Universite Paul Sabatier, 31062 Toulouse, France 



Motivated by the observation of a spin-glass transition in almost disorder-free Kagome antiferromag- 
nets, and by the specific form of the effective low-energy model of the S=l/2, trimerized Kagome 
antiferromagnet, we investigate the possibility to obtain a spin-glass behavior in two-component, 
disorder-free models. We concentrate on a toy-model, a modified Ashkin- Teller model in a magnetic 
field that couples only to one species of spins, for which we prove that a dynamic spin-glass behavior 
occurs. The dynamics of the magnetization is closely related to that of the underlying Ising model 
in zero field in which spins and pseudo-spins are intimately coupled. The spin-glass like history 
dependence of the magnetization is a consequence of the ageing of the underlying Ising model. 



PACS Nos : 75.10.Jm, 75.10.Nr, 75.50.Lk 



I. INTRODUCTION 

The experimental investigation of quantum magnets 
is undergoing rapid progress with the synthesis of new 
and better controlled samples. Of special interest are 
frustrated magnets, for which very unusual behavior has 
been reported Q|. In particular, a behavior reminiscent 
of spin glasses has been reported in a number of Kagome 
antiferromagnets ( AF) jl],|| . While the role of disorder in 
these phenomena cannot be excluded yet, the persistence 
of this behavior in progressively cleaner samples calls for 
explanations in terms of disorder- free models. 

The possibility of a spin-glass behavior in the spin 1/2 
Heisenberg model on the Kagome lattice without any 
disorder was first suggested in the pioneering work of 
Chandra et al ^ in 1993. Since the classical Heisenberg 
model on the Kagome lattice does not exhibit a spin- 
glass behavior, the main problem one is facing is how 
to include the quantum aspects of the problem into a 
description in terms of classical variables for which one 
can use standard techniques to study the dynamics. In 
Ref. H , the authors assumed that quantum fluctuations 
select coplanar configurations, which led them to con- 
centrate on the anisotropic XY version of the model. For 
this model, they suggested the presence at finite temper- 
ature of 3-spin order, and of a glassy transition related to 
the binding of non-abelian defects that are expected to 
be the natural point defects associated with this type of 
order. For details, the reader should consult Ref. j|. To 
prove or disprove this scenario turned out to be difficult, 
and whether this indeed provides an explanation of the 
spin-glass behavior of Kagome antiferromagnets without 
introducing any disorder remains unsettled. Since then, 
clear evidence of glassiness has been reported by Chandra 
et al 101 for another class of disorder-free models describ- 



ing periodic Josephson arrays in a transverse magnetic 
field, but these models are not directly related to the 
Kagome AF. 

In parallel, a lot of progress has been made in the un- 
derstanding of the spectrum of the S=l/2 AF Heisen- 
berg model on various frustrated lattices, in particular 
the Kagome and the pyrochlore ones |t). So far, it is 
well established that frustration can have two effects: It 
can open a gap to triplet excitations |)| , like for the non- 
frustrated spin 1 chain |^| , but it can also lead to a pro- 
liferation of low-lying singlets inside this gap, like for the 
S=l/2 Kagome antiferromagnet 0. These singlets can 
be interpreted as RVB (Resonating Valence Bond) states 
[p[-^0[, and they could lead to a power-law behavior of 
the low-temperature specific heat [[ll]-[l3| . Experimental 
systems known so far have a larger spin however (3/2, 
5/2,...) jjjj, and there is room for new physics in these 
cases since the presence of a singlet-triplet gap is un- 
likely given the rather small value already reported for 
S=l/2. Possible implications of this strange spectrum re- 
garding in particular a possible spin-glass behavior have 
not been discussed yet, mainly due to the lack of methods 
to attack this problem. 

In this paper, we continue the quest initiated in Ref. 
H for spin-glass behavior in the disorder-free, quan- 
tum Heisenberg model on the Kagome lattice. However, 
building on the recent results obtained on the low-energy 
spectrum of the model, we propose another approximate 
way of including quantum fluctuations into a classical 
description. The starting point is the effective model ob- 
tained in Ref. Jl4| for the spin 1 /2 Heisenberg model on 
the trimerized Kagome lattice (see Fig. (0)). This is 
a modified version of the spin 1/2 Heisenberg model on 
the Kagome lattice in which the exchange integrals take 
two different values J and J' according to the pattern of 
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Fig. ([!]). This is actually the relevant description of the 
Kagome layers in SrCrg p Gai2-9 P Oi9 since the presence 
of a triangular layer between pairs of Kagome layers lead 
to two types of bonds with precisely the pattern of Fig. 
(|l|) . Then, since the ground state of a triangle is fourfold 
degenerate and can be described by two spin 1/2 degrees 
of freedom, the total spin a and the chirality pscudospin 
t, it was shown in Ref. [Q that the low-energy effective 
Hamiltonian in the limit J' <C J can be written 

H=(j'/9) U'^U'r Ki= E ^ 

<i,3> <i,0> 

% = (1 - 2( aij rr + 4-t+))(1 - 2(/Vr + $t/)) (1) 

where < i, j > denotes pairs of nearest neighbors. In HJ-, 
ctij and Pij are complex parameters that take the values 
1, exp(27ri/3) or exp(47ri/3) depending on the bond (for 
details, see Ref. 141). 




FIG. 1. Trimerized Kagome lattice showing the difference 
between strong bonds (black triangles) and weak bonds (white 
triangles). 

In the case of the trimerized Kagome lattice, quan- 
tum fluctuations thus introduce an extra local degree of 
freedom in the Hamiltonian. For the general case, the 
derivation is no longer valid since J' /J is not a small 
parameter, but the low-energy spectrum found in ex- 
act diagonalizations of small clusters is still consistent 
with such a description: The exponentially large number 
of low-lying singlets suggest that there must be another 
quantum number in addition to the spin to classify these 
eigenstates. This is precisely the role played by the chi- 
rality pseudo-spin in the trimerized Kagome lattice. A 
detailed comparison of both spectra has actually shown 
that they are very similar indeed [ fiol . 

We would therefore like to propose a new approach 
to the spin-glass behavior of disorder-free quantum mag- 
nets in terms of models with two-local degrees of free- 
dom. Our basic motivation for concentrating on such 
models is that they are a priori good candidates to ex- 
hibit spin-glass behavior even without any disorder in the 
Hamiltonian. The basic argument is the following. In a 
field-cooled (FC) experiment, in which one quenches the 
sample from high temperature in a magnetic field, the 
system will automatically choose configurations in which 



the spins are polarized, and the pseudo-spins will adapt 
to keep the energy as low as possible. Hence the result- 
ing magnetization at low temperature is expected to be 
significant. However, in a zero-field cooled (ZFC) exper- 
iment, the system is first quenched without a magnetic 
field, and it will use both degrees of freedom to minimize 
the energy. When the magnetic field is switched on, in 
order to polarize, the system will have to overcome the 
intimate mixing of spins and pseudo-spins to get a signif- 
icant polarization. This is likely to take a very long time, 
and the apparent susceptibility will be much smaller than 
the one measured in a FC experiment. 

The analysis of the model of Eq. ([!]) with the spin 
and the pseudo-spin treated as classical Heisenberg vari- 
ables is a considerable task, and before starting such an 
endeavor, one would like to know whether the scenario 
outlined in the previous paragraph can indeed lead to a 
spin glass behavior. In fact, other models, like the fully 
frustrated XY model Jll| or some vector models ]l6| , can 
be described in terms of spin and chirality variables, and 
no spin-glass behavior was ever reported for these mod- 
els. Consequently the rest of the paper is devoted to a 
detailed analysis of a toy model to test whether the pres- 
ence of two local degrees of freedom can indeed lead to 
a spin-glass behavior. The simplest model of this kind is 
a modified Ashkin- Teller model defined by the Hamilto- 
nian: 

H = JJ2 SiSjTiTj -hY, S * (2) 

In this Hamiltonian, Si and Ti are Ising variables that 
describe the spin and the pseudo-spin respectively. Note 
that the magnetic field h is coupled only to the spin de- 
gree of freedom. Since frustration has already been in- 
cluded in the model as an extra degree of freedom, there 
is no need to work on a frustrated lattice any more, and 
for simplicity we study this model on a square lattice. 
This model is of course a very dramatic simplification 
since the Heisenberg spins and pseudo-spins are replaced 
by Ising variables. Still, as we shall see, the physics of 
this model is very rich, and to a large extent confirms the 
simple picture of the previous paragraph. Shortcomings 
that might be overcome by going to the more physical de- 
scription in terms of Heisenberg spins will be discussed 
in the last section. 



II. MONTE-CARLO RESULTS 

Let us first discuss the equilibrium properties of this 
model. If h = 0, this model is equivalent to the antifer- 
romagnetic Ising model after the local gauge transforma- 
tion o~i = SiTi, and all states have an additional degen- 
eracy of 2 N " Ue " since Oi = 1 (resp. —1) can be achieved 
with (Si, 2;-) = (1,1) and (-1,-1) (resp. (1,-1) and 
(— 1, 1)). In fact, for any value of h, the partition func- 
tion can be factorized as Z — Zi s - mg Zs where Zi s i ng is the 
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partition function of the Ising model on a square lattice, 
and Zs is the partition function of paramagnetic spins in 
a magnetic field. Accordingly, the free energy per site / 
is given by: 



/ = /i S mg-fc B Tln(2cosh(/?/ l )) 



(3) 



with /3 = 1/ksT. The magnetization per site is defined 
as usually by m = —df/dh, and since f ising does not 
depend on h, we obtain m = tanh(/3/i). So if the sys- 
tem has reached its equilibrium state, the magnetization 
smoothly saturates as the temperature goes to zero. 

However, as we shall see below, this equilibrium state 
might be very difficult to reach depending on the history 
of the system. To be specific, let us consider the follow- 
ing protocol: The quench takes place at t = 0, the field 
is switched on after a waiting time t w , and the magne- 
tization is measured after a measuring time t m elapsed 
after the field has been switched on. The time elapsed 
between the quench and the measurement is denoted by 
t. In a FC experiment, t w = and t = t m , while in a 
ZFC one, t w > and t = t w + t rn . 



Equilibrium 
m=tanh((3h) 




ZFC, t =500, t =500 



FIG. 2. Temperature dependence of the magnetization for 
different measurement conditions. 

To mimic such experiments, we have performed Monte 
Carlo simulations of the model of Eq. (|^). The elemen- 
tary step consists in flipping either a spin or a pseudo- 
spin according to Glauber prescription. The site and the 
variable (spin or pseudo spin) are chosen randomly, and 
the time unit corresponds to a number of steps equal to 
the number of degrees of freedom, i.e. twice the number 
of sites. In all numerical experiments reported below, 
the starting configuration is completely random, corre- 
sponding to an infinite temperature, and when it is not 
switched off, the magnetic field h is equal to 0.2 in units 
of J. Typical results for the magnetization as a function 
of temperature are shown in Fig. (Q) for a system of 
size 400 x 400. Below a temperature T g which is smaller 
than the Curie temperature T c ~ 2.27 J of the underlying 
Ising model, there is a clear difference between the ZFC 
and FC measurements: The ZFC magnetization drops 
quite abruptly, as in typical spin-glasses. Note that these 



curves depend on t m and t w (this dependence should of 
course vanish if t m was infinite) but only weakly in a 
large parameter range, as we shall explain below. 

To gain some insight into the origin of this behavior, 
it is very useful to study the time evolution of the mag- 
netization at fixed temperature for different values of t w . 
Typical results are given in Fig. (|J). The most salient 
features are: i) A much faster increase at short times for 
the FC experiment (t w = 0) than for the ZFC ones; ii) A 
similarity of the shape of the ZFC curves. In fact, a very 
good scaling can be obtained if we plot the magnetiza- 
tion as a function of t m /t w (see inset of Fig. (||)). Such 
a scaling is typical of ageing phenomena |17j . 




FIG. 3. Time dependence of the magnetization at given 
temperature for various measurement conditions (FC and 
ZFC with t w = 500, 1000, 2000, 3000). Inset: Plot of the ZFC 
data as a function of tm/ tw- 



ill. ANALYTICAL APPROACH: AGEING AND 
PERSISTENCE 

To understand the behavior of the ZFC and FC experi- 
ments let us first consider the dynamics of the underlying 
Ising model in the absence of a magnetic field. After a 
quench below the Curie temperature T c , the model has 
a spontaneous staggered magnetization in the variables 
Oi in terms of statics. However, in the absence of a field 
that breaks the symmetry between the two ground states, 
no global magnetization develops for an infinite system 
and the system stays out of equilibrium. Domains of the 
two phases form and coarsen with a characteristic length 
scale y/t jl8| as the evolution is via diffusion of the do- 
main walls and coalescence of domains. The ergodic time 
t eTg for such a system scales with the linear size L of the 
system like L 2 , and thus as long as t <C t CIg , the system is 
out of equilibrium and has no global magnetization. We 
have checked that for T = .6J and L = 400, the ergodic 
time is of the order of 4000. Since this regime is the only 
one accessible for very large samples, we will limit our 
discussion to that regime. 
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Another point worth noticing about the dynamics of 
the Ising model is that, after the short initial stage where 
spins with 3 or 4 parallel neighbors have disappeared, 
flips take place only at sites with two neighbors up and 
two neighbors down. For such sites, the effective field 
coming from the coupling to the neighbors is equal to 
zero, and we call them 0-sites. Since there is no energy 
cost to flip the spin Cj of such a site, the probability to 
flip according to Glauber dynamics is 1/2. So the flipping 
rate should be half the concentration of 0-sites. We have 
checked numerically that this is indeed the case after a 
few sweeps. For T = .6J, this is true as soon as t > 40. 

Let us now turn to an analysis of the problem in the 
presence of a magnetic field. Let us forget the first few 
steps and concentrate on times where flips take place at 
0-sites. If the magnetic field is small compared to J, a 
condition usually fulfilled in experiments, the dynamics 
of the Ising spins o~i will be essentially unaffected, and 
the magnetic field will just favor configurations where 
the spins Si are parallel to the field Now let us 

suppose that a 0-site remains so for a while. The equation 
that governs the appearance of a magnetization for the 
S spins can be easily deduced. If we denote by p+ (resp. 
p_) the probability for S to be up (resp. down), Glauber 
dynamics implies that p+ satisfies the equation 

dp+ _ e 0h e~ ph 

-fa ~ e -Hh + e ph P~ ~ e -ph + e ph P+ W 

The magnetization, which is related to p+ and p_ by 
p + = (1 + m)/2 and p_ = (1 — m)/2, is thus given by 

m = tanh(/3/i)(l - exp(-t/r)) (5) 

where the relaxation time r is equal to 1 in the chosen 
time units. This is a very short time, especially consid- 
ering the very small flipping rate which is achieved af- 
ter a few steps (already below 10 -2 after 100 sweeps for 
T = .6 J). Under these circumstances, the sites that have 
been 0-sites in the presence of the magnetic field should 
have enough time to reach equilibrium, and on average 
their magnetization should be equal to tanh(/3/i). So if 
we call Co(t m ,t w ) the proportion of sites that have been 
0-sites between t w and t w +t m this simple argument leads 
to the prediction that 

m(t m ,t w ) =tanh(Ph)co(t m ,t w ) (6) 

Note that co(t m , t w ) depends only on the dynamics of the 
underlying Ising model and not on the magnetic field. 
To check this prediction, we have calculated co(t rn ,t w ) 
for different t w corresponding to our ZFC numerical ex- 
periments. The agreement is very good - the curves are 
indistinguishable from the the ZFC calculations of m on 
the scale of Fig. (||) - and we have checked that it remains 
so as long as t w is not too small so that the magnetization 
is indeed controlled by 0-sites. This analysis shows that 
the ageing behavior of the model of Eq. (1) is indeed 
closely related to the dynamics of the underlying Ising 
model. 



To understand the magnetization process we therefore 
only have to consider Co(t m , t w ) for the Ising model. After 
the first few steps, domains can be identified, and 0-sites 
are on the boundary of the domains. On a square lattice, 
they correspond to the diagonal portions of the bound- 
ary. Now, after a time t w the characteristic size of the do- 
mains is of order y/t^. Consequently the total number of 
domains is of order 1 jt w , and the total length of domain 
walls is proportional to l/s/t^. When t m <C t w we re- 
mark two important points: (i) the domain walls appear 
locally flat, the effects of surface tension and hence cur- 
vature are negligible and hence the domain walls diffuse 
(as zero modes are the only ones generating dynamics by 
this stage), a given point on a domain wall will therefore 
diffuse a length \/*m; (ii) one can neglect the interaction 
between domain walls, that is to say the coalescence of 
domains. The total number of spins which flip between 
t m and t m + t w is therefore proportional to v / t^x total 
length of domain walls. We therefore find that at short 
times t m (compared to t w ) one has 

ca(t mi t w ) = Const. J (7) 

V "w 

where the constant should be independent of t w . This 
provides an excellent fit of the data at short times (see 
Fig. (|). We have checked that Const. ~ 0.05 at T = .6 J 
is indeed independent of t w . This argument provides a 
very simple explanation of the fact that the larger t w , the 
slower the initial increase in the magnetization. 
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FIG. 4. Fits of the FC and ZFC data at T/J = 0.6 with 
the help of Eqs. (7-9). 

In the case where t w = 0, the above argument can- 
not be applied as it is because of the initial steps, where 
flips occur at sites with 3 and 4 parallel neighbors as well, 
and because domains are not present right away. Still the 
same kind of reasoning suggest that every spin which has 
flipped has an average magnetization tanh(/3/i). There- 
fore one has m(t, 0) = tanh(/3/i)(l — p(t, 0)) where p(t,0) 
is the probability that a given spin has not flipped before 
time t. The quantity p(t ) 0) has received much atten- 
tion in the literature [£2| and is called the probability of 
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persistence. Extensive numerical simulations have shown 
that for large times pit) ~ 1 /t 6 , where 6 is the persistence 
exponent which has been measured as ~ 0.2 |23] for the 
two dimensional Ising model. One therefore deduces that 

m(t,0) = tanh(/3/i)(l - (r/t) e ) (8) 

Here r is a microscopic time scale related to the flipping 
rate of the spins. This provides and excellent fit of the 
data except at very short times (see Fig. [|). 

In the case of t w nonzero, and for t t w , one may 
adopt a coarse grained point of view and cut the sys- 
tem up into units of size \ft w . These blocks may be re- 
garded as effective spins of the sign of the average block 
staggered magnetization. This is related to block per- 
sistence introduced in pl| . Moreover, the characteristic 
time scale r(t w ) for these effective spins is proportional 
to the time it takes to reverse the staggered magnetiza- 
tion of a given block, since this is done by diffusion one 
has r(t w ) ~ t w . One may now apply the reasoning of the 
case t w = to obtain 

m(t, t w ) = tanh(/3/i)(l - (r(t w )/t) e ) (9) 

for t 3> t w . This again provides a reasonable fit of the 
data (see Fig. (^)), and we have checked that r(t w ) is 
proportional to t w . 

To summarize, the very slow onset of magnetization 
in ZFC experiments follows scaling laws typical of age- 
ing, and this is the origin of the spin-glass behavior. 
This ageing is a consequence of the very slow march to- 
wards equilibrium of the underlying model in which spins 
and pseudo-spins are coupled into an effective Ising spins 
since the physical field, which couples only to the spin 
and not to the pseudo-spin, does not act as a symmetry 
breaking field for that model. 

IV. DISCUSSION 

While the difference between FC and ZFC measure- 
ments of the magnetization is very typical of spin-glasses, 
a bona fide spin-glass has many other characteristics 
which are not shared by the present model. 

First of all, for an infinite system the freezing tran- 
sition of a standard spin glass is believed to be a ther- 
modynamic transition, and the freezing temperature is 
expected to be independent of the protocol. In the 
present case, in order to flip, a spin inside a domain must 
be able to cross an energy barrier of 8 J. The Arrhe- 
nius law gives therefore a characteristic flipping time of 
r Q ~ exp(8J/T). Such activated flips are only observed 
after a time of measurement of order r a . One may there- 
fore define a temperature depending on the time scale 
of the measurement T g {t m ), such that if r a ^> t m such 
activated spin flips do not occur and hence the equilibra- 
tion of the magnetization within domains is not possible. 
The crossover between the two regimes is therefore given 



by T g (t m )/J r-j 8/ln(t m ). So, unlike the standard ther- 
modynamic spin glass transition, the freezing tempera- 
ture depends on t m . This estimate of T g agrees with our 
Monte Carlo results. If one takes a measurement time of 
t m = 500, as for the curves shown in Fig. (|J), one finds 
that T g / J ~ 1.287, close to the value of T g shown in Fig. 

(§■ 

Another original aspect of our data with respect to 
standard spin glasses is the minimum of the ZFC mag- 
netization at a temperature far below T g . To under- 
stand this, one may make an approximate decomposi- 
tion of the magnetization in terms of the activated com- 
ponent m a within the domains and the component gen- 
erated by the domain walls m c i- We find m a (t m ,t w ) — 
tanh(/3/i)(l — exp(— a(h)t m /T a )) (solving the Glauber dy- 
namics explicitly for a single spin with its four neighbors 
fixed parallel shows that a = 1 + cosh(2/3/i) + 0(l/r„)). 
As m a is generated within domains and hence in the bulk 
of the system (that is to say not at domain interfaces), 
it is very weakly dependent on t w . Another type of ac- 
tivated process may occur at domain interfaces where 
a spin with three neighbors antiparallel and one parallel 
flips. The energy barrier for this is 4 J. This gives another 
characteristic time r* with a corresponding time depen- 
dent temperature T*(t m )/J ~ 4/ln(t m ). For t m = 500 
as in Fig. (Q) one finds T*/J — 0.6436 which corresponds 
to the minimum seen on the ZFC curve. For T -C T*, 
the system behaves as if J is infinite on the experimen- 
tal time scale and only domain wall diffusion occurs, the 
only energy involved being the external field energy h. 
Here the dynamics is well described by Eq. (||), hence 
the measured magnetization m(t m ,t w ) increases as T is 
lowered below T*. 

Both these effects are dynamic in nature and show that 
the behavior we have observed is typical of the out of 
equilibrium dynamics observed in spin-glasses. This con- 
clusion is in fact consistent with another aspect of our 
data, namely the sensitivity of the results to the dynam- 
ics used in the simulation. To perform Monte Carlo sim- 
ulations, we have made the assumption that spins and 
pseudo-spins flip independently of each other. However 
if one allows simultaneous flips as well, the behavior will 
be very different: These moves will allow the magnetiza- 
tion to develop inside the Ising domains, and the system 
will polarize after only a few sweeps. In other words, the 
breaking of ergodicity is related to the dynamics. 

Another way to check how different the present model 
is from standard spin glasses is to study the tempera- 
ture dependence of the non-linear susceptibility, which is 
expected to diverge at the transition in a spin-glass. Pre- 
liminary results p3] indicate that the non-linear suscep- 
tibility is not singular in the present case. The absence 
of another thermodynamic singularity below the Curie 
temperature T c is also clear from Eq. (||). 

However, other characteristic aspects of spin-glasses 
are also shared by the present model. For instance, 
the thermo-remnant magnetization tutmr shows a very 
strong dependence on the time t w elapsed between the 
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quench and the switching off of the magnetic field: The 
longer t w , the slower the decay. A detailed analysis of 
these results along the same line is in progress. 

Coming back to our original purpose, namely the ex- 
planation of the low-temperature behavior of Kagome an- 
tiferromagnets, the present model has both merits and 
drawbacks. The very clear difference between FC and 
ZFC magnetizations is the most interesting aspect of the 
results. It shows that the presence of an extra degree 
of freedom which is not coupled to the external field 
can indeed lead to a glassy behavior at low tempera- 
tures. However most of the differences with standard 
spin- glasses are problematic in that respect: No experi- 
mental indication of a dependence of T g on the protocol 
was reported, and the non-linear susceptibility is indeed 
enhanced close to T g . This is not a final blow however. 
In fact, these differences with standard spin glasses all 
depend on the fact that flipping simultaneously a spin 
and a pseudo-spin leaves the Ising spin unchanged. This 
symmetry will not be present in more realistic models 
where spins and pseudo-spins are treated as Heisenberg 
variables, while the underlying mechanism for the diffi- 
culty that the system will have to magnetize after be- 
ing cooled in zero-field is still expected to apply. More 
realistic models with two degrees of freedom treated as 
Heisenberg spins are therefore good candidates for effec- 
tive models to get a spin-glass behavior in disorder free 
magnets. Work is in progress along these lines. 
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